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We show that the two-band nature of noncentrosymmetric superconductors leads naturally to 
an anomalous Josephson current appearing at zero phase difference in a clean noncentrosymmetric 
superconductor/ferromagnet/noncentrosymmetric superconductor junction. The two-band nature 
provides two sets of Andreev bound states which carry two supercurrents with different amplitudes. 

When the magnetization direction of the ferromagnet is suitably chosen, two supercurrents experi¬ 
ence opposite phase shifts from the conventional sinusoidal current-phase relation. Then the total 
Josephson current results in a continuously tunable ground-state phase difference by adjusting the 
ferromagnet parameters and the triplet-singlet ratio of noncentrosymmetric superconductors. The 
physics picture and analytical results are given on the basis of the s+p wave, while the numerical 
results are reported on both s-fp and d+p waves. For the d+p wave, we find novel states in which 
the supercurrents are totally carried by continuous propagating states instead of discrete Andreev 
bound states. Instead of carrying supercurrent, the Andreev bound states which here only appear 
above the Fermi energy block the supercurrent flowing along the opposite direction. These novel 
states advance the understaning of the relation between Andreev bound states and the Josephson 
current. And the ground-state phase difference serves as a tool to determine the triplet-singlet ratio 
of noncentrosymmetric superconductors. 

PACS numbers: 74.50.-|-r, 74.70.Tx, 74.20.Rp 


I. INTRODUCTION 

The noncentrosymmetric superconductor (NCS) has 
attracted much attention for the coexistence of spin- 
singlet and spin-triplet superconductivity and the pos¬ 
sibility of topologically nontrivial surface statesAZ. The 
unusual properties of NCS originate from the absence 
of inversion symmetry from the crystal structure, which 
permits an antisymmetric spin-orbit coupling (SOC) odd 
in electron momentum, and leads to a chiral ground 
state. The SOC mixes the spin-singlet (even-parity) com¬ 
ponent and the spin-triplet (odd-parity) component in 
the superconducting pairing potentials^— . The list of 
NCSs has grown to include dozens of materials such as 
Li 2 Pda;Pt 3 _a,BSiii 2 , Y 2 C 3 SS, and the heavy-fermion com¬ 
pounds CePtsSiiS, CeRhSigS, and CelrSigSS. These com¬ 
pounds have various crystal structures and hence various 
forms of SOC and mixed triplet-singlet pair symmetry. 
The actual superconducting pairing mechanics and pair¬ 
ing state symmetry realized in these NCSs are still un¬ 
clear. 

In CePtsSi and several other Ce-based NCSs, the spin- 
singlet and the spin-triplet components are expected to 
appear in comparable magnitudes'^, which may lead to 
more exotic effects. As a result, the NCS has two effec¬ 
tive superconducting gaps that are relevant to two sets of 
spin-split Andreev bound states. The relative magnitude 
of the two parity components determines the relative size 
of the two gaps and the main property of a NCS. The 
question of how to determine the triplet-singlet ratio has 
attracted many theoretical and experimental efforts^^^— . 
As one of the main experimental methods in the field of 
superconductivity, the Josephson effect has been widely 


studied in junctions based on NCSs to investigate the 
characteristic triplet-singlet ratios. There were efforts 
focusing on the investigations of the steps in the current- 
voltage characteristics^^, the low-temperature anomaly 
in the critical current^^, and the transition from a 0 
junction to a 7r/2 junctioiJ^ in Josephson junctions with 
NCSs. The magnetic Josephson junction with NCSs has 
also been studiedi^ and the authors focused on the high- 
order harmonics in the charge and spin current-phase re¬ 
lations and the possibility of O-tt transitions. 

On the other hand, much attention has also been paid 
to the anomalous Josephson effect which means a ipQ- 
junction with arbitrary ground-state phase difference (po 
other than 0 or tt. Usually, the supercurrent in a Joseph¬ 
son junction vanishes when the phase difference between 
the two superconductors becomes zero and the current- 
phase relation (CPR) is sinusoidal I{ip) = /cSin(<p) in 
the tunnelling limit^^. Whereas an anomalous Joseph¬ 
son current la flowing even at zero phase difference 
{(fi = 0) has recently been predicted in various types of 
Josephson junctions^i^i^. The anomalous supercurrent 
is equal to a phase shift (po in the conventional CPR, 
i.e., I{(p) = IcSm((p — (po). In general, there are two 
prerequisites for the emergence of a tpo-junction: (i) two 
sets of spin-split Andreev bound states (ABS) with op¬ 
posite phase shifts iipo compared with the conventional 
CPR, and (ii) different amplitudes of the supercurrents 
carried by the two sets of ABS. Especially, the Josephson 
junction formed on the surface of topological insulators 
by the proximity effect can be considered as the limiting 
case where only one set of ABS remain^^. 

In Josephson junctions based on NCSs, the second 
prerequisite is naturally met due to the two-band 
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FIG. 1: (Color online) Schematic diagram of the NCS/F/NCS 
junction. The ferromagnetic layer F has a finite width L. The 
exchange field h lies in the x-y plane and makes an angle 6 
with the a:-axis. 


nature of superconductivity in NCSs. The ABS are 
naturally spin-split into two sets and the two gaps with 
different sizes ensure that the supercurrents carried 
by two sets of ABS have different amplitudes. To 
meet the first prerequisite, we introduce a ferromagnet 
into the Josephson junction which brings phase shifts 
to ABS. In this study, we investigate the anomalous 
Josephson effect in a noncentrosymmetric superconduc- 
tor/ferromagnet/noncentrosymmetric (NCS/F/NCS) 
junction. It is shown that the ground-state phase differ¬ 
ence 1^0 is sensitive to the triplet-singlet ratio of NCSs. 
Therefore, the anomalous Josephson effect serves as a 
mechanism to determine the unknown triplet-singlet 
ratio of a NCS. The physics picture and analytical 
results are given on the basis of the s+p wave, while 
the numerical results are reported on both s+p and d+p 
waves. For d+p wave, we find novel states in which 
the supercurrents are totally carried by continuous 
propagating states instead of discrete Andreev bound 
states. 

The paper is organized as follows. In Sec. II we 
present the model Hamiltonian and introduce the nu¬ 
merical method based on the lattice Green’s function to 
solve the CPR. In Sec. Ill we present the analytical re¬ 
sults of the normal incidence component in the case of 
s-bp wave. The numerical results and relevant discussion 
on three types of pair potentials s-bp, d-bp, d-bf will be 
given in Sec. IV. Finally, the conclusion will be given in 
Sec. IV. 


II. MODEL AND NUMERICAL METHODS 


We consider a two-dimensional NCS/F/NCS junction 
in the clean limit. A schematic diagram of the junction 
under study is shown in Fig. [I] The ferromagnetic layer 
F has a finite width L, and an exchange field h whose 
direction is in the x-y plane and makes an angle 9 with 
the x-axis. 

The numerical method used to evaluate the supercur¬ 
rent is the lattice Green’s function technique. We con¬ 
sider the junction in a square lattice with the lattice con¬ 
stant a. The lattice lies in the x-y plane. The ferro¬ 
magnetic layer F sandwiched by the two NCS electrodes 
locates in the region [0, L = Na] with N is the number 
of columns. During the tunnelling processes in the x- 
direction, the transverse momentum ky is assumed to be 
conserved. In this context, the Hamiltonian of the hybrid 
junction reads 


^ ~ ^ [^kky d] cL,fc„G(T,/c, 

1(7,ky 


1(7,ky 


+ jhlapa' + Xlapcr') cr' ,ky 

1,1' ,ky (7(7' 

1,1',ky (7(7' 


( 1 ) 


where ^ (Q(t,/c„) is the creation (annihilation) opera¬ 
tor of an electron in column I with spin cr (= f or j,) and 
transverse momentum ky. The on-site energy ei^ky has 
the form —2t cos Ky with Ky = kya. The Fermi energy p, 
and the nearest-neighbor hopping integral t are assumed 
to be the same in the whole junction. The hopping coeffi¬ 
cients of the exchange interaction and the Rashba 

SOC are given by the matrices in the spin space 


hw = h • crSip 

\u’ = \{cFxS\nKy5ip+ iay5i^ip/2) ( 2 ) 

where h and A denote the exchange field and the Rashba 
strength respectively. The pair potential coefficient 
readsiS^ 


J 




A [(g -b ysinKyax) Sip + ^qaySi^ip] iay, 

A {[±§9 sin ATy ± iq{l — cosKy)ax + qs\nKyay\ Si^ip — 2qsv[\Kyay5ip) iay, 

A [{2cosKy5ip - Si^ip) {q + qsinKyax) + q{+i cos KySi^^zip ± ^Si^ 2 p)cky\ i(Xy, 


s + p-wave 
dxy + p-wave 
dx 2 -y 2 -bp-wave 


(3) 


where A is the BCS gap function which takes Aq at zero simplicity, the pair potentials in the two NCS leads are 
temperature and vanishes at critical temperature Tc. For set to be equal in amplitude, ai with i = x, y, z are the 













3 


Pauli matrices. The parameter q {q = 1 —(?) characterizes 
the percentage of spin-singlet (triplet) component which 
reads qL (qh) and qn (qji) in the left and right NCS re¬ 
spectively. In Eq. we also omit the phase difference 
between the two NCSs for simplicity. The phase differ¬ 
ence is set to be tp. Note that the exchange field exists 
only in the F layer while the Rashba SOC and the pair 
potential only exist in the two NCSs. 

In the F region (0 <C ? IV), the charge operator in 
column I with momentum ky is defined as 

Pi,ky = ec\j,^aQCi,ky, ( 4 ) 

where ci,ky = [cit,fey > Qt.fey ^ o'o is the unite ma¬ 
trix, and t is the time. By using the Heisenberger equa¬ 
tion ihdfp = [p, H], the operator of supercurrent is found. 
Then we can construct the Green’s function to calculate 
the supercurrent I through column I as follows^ 

/ = ^ y Tr [tieG< (/, I - 1) - etG< {I - 1, 1)] dKy (5) 

where i = — trs (g) crp and e = —era ® uo denote the 
hopping matrix and the charge matrix respectively. T 3 is 
the Pauli matrix in Nambu space and e (> 0) is the unit 
charge. In equilibrium, the lesser-than Green’s function 
G^ equals 


G< = y'^/(F)[G“-Gn (6) 

where / {E) is the Fermi-Dirac distribution function. 
The retarded (advanced) Green’s function (G“) can 
be numerically calculated by the recursive method. 

Besides the supercurrent, the ABS spectra can be also 
calculated numerically. It is known that the ABS results 
in the peaks of particle density within the superconduct¬ 
ing gap. Therefore, by searching the peaks of particle 
density in column I {N ^ Z ^ 1) 

pi = --Ini[Tr{G^l,l)}] (7) 

TT 

at a given phase difference ip, the energies of ABS can 
be located. Then we scan ip and obtain the ABS spec¬ 
trum which is useful for understanding the behavior of 
supercurrent. 


III. ANALYTICAL RESULTS 


Before we discuss our numerical results, we present 
the analytical results of the normal incidence component 
in the s+p wave case, which is very helpful to under¬ 
stand how the junction under consideration becomes a 
(/ 3 o-junction. We start with the Bogoliubov-de Gennes 
(BdG) Hamiltonian for a NCS in the momentum space 


f eic + Aik ■ cr A(k) 

V At(k) -ek + Alk-a* 


( 8 ) 


TT fd{kl+kl) t2t2 , 

Here, Ek = --^ ^ - M is the spin- 

independent part of band dispersion with p the chemical 
potential, and a is the vector of Pauli matrices. Aik • cr 
is the antisymmetric SOC with Ik = kyX — kxY and A 
the SOC strength. The superconducting gap function is 
A(k) = /(k)(As -I- Atdk • (j)i(jy, where As = Aq? and 
At = Ao(l — 9 ) are spin-singlet and spin-triplet super¬ 
conducting gaps respectively with Ag = Ag -I- At and q 
turns between purely spin-triplet (q = 0) and purely spin- 
singlet {q — 1) pairings. We assume Ag, At, and Ag are 
positive constants, and the orbital-angular-momentum 
pairing state is described by the structure factor /(k). 
In the case of s-l-p wave, /(k) = I. The spin-triplet pair¬ 
ing vector is aligned with the polarization vector of the 
SOC dk = Ik/^F where kp is the Fermi wave vector and 
taken as the unit of the wave vector. When the SOC 
splitting is much less than the chemical potential, we can 
use the Andreev approximation kp kp = 1 with kp 
the spin-split Fermi wave vectors. 

To diagonalize the kinetic term, it is convenient to ex¬ 
press the Hamiltonian in the so-called helicity basis. We 
introduce the following spin rotation transformation 


f U 0 

Vo 



1 

1 ) ’ 


(9) 


where (j) = arctan^ky/kx) is the incident angle of quasi¬ 
particles. Under the rotation R, the normal part of the 
Hamiltonian is diagonalized and the superconducting gap 
function is transformed to 


C/A(k) ([/•)-■ (10) 


with /(k) = —ie“*‘^/(k). The factor is the same 

for the left and right NCSs and thus has no net effect on 
the Josephson effect in the first harmonic approximation 
where the normal reflection is absent. The Hamiltonian 
of a NCS in the hilicity basis reads 


/ Ek — Afc 0 0 

0 Ek + A/c A+/(k) 
0 A+/(k) -Ek - \k 

\ A_/(k) 0 0 


A-/(k) \ 
0 
0 

—Ek T \k j 


( 11 ) 


where A± = At ± Ag. It is clear that the Hamiltonian 
shows a two-band nature that there are two bands with 
different superconducting gaps |A±| are uncoupled in the 
helicity basis. One band is for the Cooper pair made 
of spin-up electron and spin-down hole with gap |A_|, 
the other band is for the pair of spin-down electron and 
spin-up hole with gap A_|_ with respect to the helicity 
basis. Since the critical Josephson current is linear in 
the gap, the two bands provide two supercurrents with 
different amplitudes. Thus the second prerequisite for 
the emergence of a (/?g-junction is naturally reached. 

The first prerequisite for a (/ 9 g-junction is easy to meet 
by utilizing a ferromagnetic interlayer. The middle F 
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layer can bring opposite phase shifts to the two supercur¬ 
rents provided by the two bands because the spin-triplet 
component of the gap function is also spin-opposite pair¬ 
ing in the helicity basis. For simplicity, we consider only 
the normal incidence component with ^ = 0 which is 
dominant in the Josephson current for the case of orbital 
s+p wave. We consider a F interlayer with the exchange 
filed strength h and the width L. The magnetization 
direction is chosen to be aligned with the polarization di¬ 
rection of the SOC, i.e., the y-direction. In this situation, 
the two helical bands keep uncoupled from each other and 
the Hamiltonian of the whole NCS/F/NCS junction for 
the two bands is respectively 


/ Sk - crXk + ah{x) A^(a;)/(k)^ \ 

V A-^{x)f{k) -Ek + crXk + ah{x) ) 


where the helicity index ct = 1 is for the pair of spin- 
up electron and spin-down hole while a = —1 is for 
the pair of spin-down electron and spin-up hole with 
respect to the helicity basis, h{x) = sgn{kx)h{x) with 
h{x) = h[Q{x) — <d{x — L)], = Arp when a = ±1 

and Aa{x) = Azp[0(—-I- 0(a; — with 

(fi the macroscopic phase difference of the two NCSs. 
The F layer is expected to bring phase shifts arj to the 
set of ABS with index a and the corresponding super¬ 
current. In the dimensionless units, the phase shift is 
approximately^ rj = {k- — k+)L k, hL with the wave vec¬ 
tors for spin-up or spin-down k± = y/W^Ti Ri k ^ h/2. 
Note that the right-going and left-going Cooper pairs 
with the same helicity index a actually have opposite 
real spin with respect to the y-direction. It is interesting 
that these two pairs have the same phase shift ar] induced 
by the F layer—, and correspond to the same gap A^ at 
the same time. 

After the two prerequisites are reached, we come to 
discuss the supercurrent carried by the two bands in the 
NCS/F/NCS junction by assuming that the ratio of the 
supercurrent to the gap is a constant for the two bands. 
In the first harmonic approximation, we have two super¬ 
currents 


with 


sin ipo 


cos (fio 


—A 2 siny 


Af cos^ y -|- A| sin^ ; 

Ai COST] 

Af cos^ T] + A^ sin^ ; 


(15) 


It is seen that the anomalous ground-state phase differ¬ 
ence ipo depends not only on the F layer induced phase 
shift r] = hL but also the triplet-singlet ratio of the two 
NCSs. Thus the arbitrary (^o-junction can be obtained 
by tuning the ferromagnet parameters and the triplet- 
singlet ratio. We can also determine the triplet-singlet 
ratio of NCSs by detecting the ground-state phase differ¬ 
ence of such a NCS/F/NCS junction. 

In the above discussion, we assume that the left and 
right NCSs have the same singlet percentage q. If the two 
NCSs have opposite q parameters, i.e., qL = ^ — qn with 
Ql{r) the q parameter of left (right) NCS, the situation 
is a little different. For the A_ band, the exchange of 
As and At changes the sign of the gap function as shown 
in Eq. CD- Therefore, there will be equivalently an 
additional phase difference tt between the left and right 
NCSs. For the A+ band, the situation keeps unchanged. 
Then the total Josephson current is 


I = —1(7 = 1 + I(7 = -l 

= P [A+ sin((/j -b 77) - |A_| sin(i 75 - 77 )] 

= 2/3 (A 2 sin(^cos? 7 -b Ai cos(/ 3 sin 77 ) 

= 2/3sin((^-(^ 0 ) (16) 


with 


sin p'q 
cos </5 q 


—Ai sin 77 


A\ cos^ 77 -b A^ sin^ ; 

A 2 COST/ 

A\ cos^ 77 -b A^ sin^; 


(17) 


For the special case of = 1 and qn = 0, the result of a 
7 r /2 junction for the triplet-ferromagnet-singlet Joseph¬ 
son junction is recovered as the same as in Ref. 


7o-=/3 I At - cr As I sin(vj - (T 77 ). (13) 

The total Josephson current is the sum I = 1^=1 +Ia=-i- 
It is noticeable that the Josephson current depends on 
only |Af — crAsI and it does not matter whether At > As 
or At < As in the first harmonic approximation. We 
refer to the bigger (smaller) one of As and At as Ai 
(A 2 ). Then the total Josephson current can be written 
as 

I = 2/3 (Ai sin 73 cos 77 -b A 2 cos 73 sin 77 ) 

= 2 / 3 sin ((/3 - (/ 3 o) (14) 


IV. NUMERICAL RESULTS AND DISCUSSION 
A. s+p wave 

For the s+p wave case, the momentum dependence of 
the gap function /(k) = 1. The magnetization direction 
of the F layer is chosen to be aligned with the //-direction 
because the contribution from normal incidence ky = 0 
is dominant. In Fig. [2l we show the numerically solved 
ABS at normal incidence in the NCS/F/NCS junction. 
With increasing ferromagnet induced phase shift 77 (from 
left to right in Fig. HD, two sets of ABS with different 
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FIG. 2: Andreev bound states at normal incidence = 0 in a NCS/F/NCS junction for s+p wave. The exchange field h in the 
F layer is aligned with the j/-direction. From left to right, the exchange field induced phase shift varies from 0 to tt with a step 
7r/4 by adjusting the exchange field strength h\ from top to bottom, q = 0,0.1, 0.3, 0.5, 0.7. Other parameters are: p — —2t, 
Ao = 0.01/i, A = O.lfi. 


helicity index a depart from each other more heavily. 
The set of ABS with index a experience a phase shift 
at] from the original degenerate ABS. With increasing 
singlet percentage q, the energy span of the set of ABS 
with helicity index <7 = 1 shrinks due to the decreas¬ 
ing minor gap Ay = A_ = At — Ag. That means the 
corresponding supercurrent with cr = 1 has a smaller am¬ 
plitude. The special case of g = 0.5 is noticeable because 
the minor gap closes completely. Only the ABS with in¬ 
dex tr = — 1 is the remaining supercurrent-carried ABS. 
Then a (/jq- junction is easily realized by taking ij = —(po, 
which agrees well with Eq. (1151) . The evolution of ABS 
with varying p and q shown in Fig. [^is fully consistent 
with the above analytical results. 

It is also noteworthy that the ABS is almost the same 
for q = 0.3 and g = 0.7 as discussed in the analytical re¬ 
sults. Because of the presence of weak normal reflections 
at NCS/F interfaces in the numerical results, there opens 
a small gap at some ABS crossing points due to the cou¬ 
pling of two ABS bands with the same spin and the oppo¬ 
site travelling direction (thus with the opposite helicity 


index a). What is interesting is that the anti-crossing 
effect for the case of At > Ag is much weaker than that 
for As > Af. We can understand this effect easily in two 
limit cases q = 0 and 9=1. For the coupled two bands 
with the opposite kx and the opposite cr, Eq. dm shows 
that At changes sign while Ag does not. That is the 
essential difference between the singlet and triplet pair¬ 
ing. Thus, the normal reflection induced coupling results 
in an anti-crossing gap opening for the triplet-dominant 
case but no obvious effect for the singlet-dominant case 
at ABS crossing points. The anti-crossing effect gives rise 
to high-order harmonics of the CPR, which is sensitive 
to the triplet-singlet ratio as reported in Ref. [l^. In this 
paper, we focus on the first harmonic of the CPR. 

The CPR of the total Josephson current is shown in 
Fig. [21 The temperature is taken as T = 0.9Tc so that 
only the first harmonic is remaining. When the F layer 
is absent, the CPR shows a normal 0-junction for all the 
values of q as shown in Fig. [3] (a). While the critical cur¬ 
rent decreases with increasing q firstly till g = 0.5 and 
then increases when q increases from 0.5 to 1. The change 
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FIG. 3: (Color online) CPR of total Josephson current for the s+p wave with various phase shift p: (a) 0, (b) 7r/4, (c) '7r/2, (d) 
37r/4. In each panel, the results for varying values of q from 0 to 1 are plotted. The temperature T = 0.9Tc. Other parameters 
are the same as those in Fig. [2] 


of the critical current coincides well with the change of 
the minor gap in Fig. because the contribution from 
normal incidence = 0 is dominant. Although the ABS 
are almost the same for q and 1 — g in the case of nor¬ 
mal incidence, the total Josephson current is not exactly 
the same for q and 1 — g because of the contribution 
from inclined incidence. For oblique incidences, the spin¬ 
splitting of the Fermi surface is enhanced and will modify 
the magnitude of triplet pairing. 

When the phase shift q = hL = 7r/4, the CPR is shown 
in Fig. 131(b) for various g. For g from 0 to 1, the ground- 
state phase difference ipo reduces firstly to the minimum 
—7r/4 at g = 0.5 and then goes back up to 0. The critical 
current also reaches its minimum at g = 0.5. All these 
features are qualitatively consistent with the ABS in Fig. 
[2]andEq. (fTSl) . Unfortunately, we cannot distinguish the 
triplet-dominant pairing from the singlet-dominant pair¬ 
ing just by the ground-state phase difference. However, 
such is not the case when the phase shift ry = 7r/2. It 
is shown in Fig. [3] (c) that the ground-state phase dif¬ 
ference of the CPR decreases monotonically from 0 to 
—TT when g goes up from 0 to I. It is important that we 
can determine the triplet-singlet ratio just by the ground- 
state phase difference of the CPR. This special feature of 
ry = 7r/2 can be understood by considering three limit 
cases g = 0,0.5,1. The case of g = 0.5 is simple. The 
gap with (7 = 1 closes and the remaining set of ABS with 


a = —1 experience a phase shift aq = —tx 12. Obviously 
we obtain a —7r/2 junction. For the cases of pure triplet 
or pure singlet pairing (g = 0, or 1), the supercurrents 
carried by two sets of ABS have the opposite phase shift 
±7r/2 and the same amplitude. Then the two supercur¬ 
rents cannel each other out. Thus the contribution to 
supercurrent from large incidence angles is dominant in¬ 
stead of from normal incidence. For large incidence angle, 
the magnetization direction (yy-direction) is not aligned 
with the spin direction of helicity basis any more. The 
singlet pairing is always opposite-spin pairing (indepen¬ 
dent of the direction of the spin quantization axis) while 
the triplet pairing is equal-spin pairing in the spin quan¬ 
tization axes perpendicular to the helicity basis. It is 
well-known that the F layer cannot bring a phase shift 
for the equal-spin pairing. And when ky ^ 0, the wave 
vector difference kx- — kx+ > /i for g = 1. So the larger 
the incident angle becomes, the less (more) the absolute 
value of the phase shift is than 7r/2 for g = 0 (1). For 
pure pairings (g = 0, or 1), the CPR should be either 
a 0-junction or a 7r-junction because of the presence of 
two supercurrents with the opposite phase shift and the 
same amplitude. Therefore the CPR tends to become a 
0-junction for g = 0 while a 7r-junction for g = 1. When 
g goes up from 0 to 1, the CPR naturally exhibits a 0- 
TT transition with smoothly changed ground-state phase 
difference. When q = 37r/4, the ground-state phase dif- 
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FIG. 4: (Color online) Particle density in the F layer at incidence angle 0 = 7r/4 for d+p wave. The singlet percentage is fixed 
to g = 0.56 to close the minor gap A_. From left to right, the exchange field induced phase shift r) = Q, 7r/4, 7r/2, 37r/4, tt; 
from top to bottom, the angle of the magnetization direction 6^ = 0, 7r/4, 7r/2, 37r/4, respectively. Other parameters are the 
same as those in Fig. O 


ference varies between —tt and —37r/4, which is similar 
to the case of r; = 7r/4 and also consistent with the ABS 
in Fig. [5] and Eq. (1X51) . The difference is that the critical 
current for q > 0.5 is now larger than that for g < 0.5 
because of the contribution from inclined incidence. 


B. d+p wave 

For the d+p wave NCS, the momentum dependence of 
the gap function /(k) = kxky/kp. The dominant contri¬ 
butions come from the two components with the incident 
angle (p = ±7r/4. We cannot choose a single magnetiza¬ 
tion direction which is aligned with both spin polariza¬ 
tion directions of these two components. That is to say, 
the coupling of two helical bands is unavoidable for at 
least one of the two components. For simplicity, we be¬ 
gan with the special case where only one gap survives and 
the other is closed. The numerical results show that this 
can occur at g = 0.56 when p = ±7r/4. The singlet per¬ 
centage q deviates from 0.5 due to the lattice model and 
the spin-splitting of Fermi surface. In Fig. 21 the particle 
density in the F layer at incidence angle </> = 7r/4 shows 
the evolution of ABS with varied strength and direction 
of the exchange field. Because the minor gap is closing, 
there are only two ABS remaining. The right-going (left¬ 
going) ABS consists of a right-going (left-going) electron 


and a left-going (right-going) hole. It is noticeable that 
two ABS have different eigen spin direction because ky is 
fixed and finite. 

When 0 = 0, the magnetization direction makes an 
angle of 7r/4 with both spin-down directions (cr = — 1 
with respect to the helical basis) of two ABS. On the one 
hand, the F layer precesses the spin of electron and hole 
by an angle p ~ hL/ cost/). Some spin-down particles 
are flipped to spin-up and enter into the NCSs without 
Andreev reflection because the spin-up gap is closed. The 
remaining spin-down particles proceed to finish the cycle 
of two Andreev reflections to form a ABS. So the particle 
density of two ABS shrinks with increasing rj. On the 
other hand, the F layer also brings a phase shift to the 
two ABS. The phase shift is opposite for right-going and 
left-going ABS. This is opposite to the previous case of 
the s+p wave where the phase shift is the same for right¬ 
going and left-going ABS. In that case, the magnetization 
direction along the y-axis is parallel to the spin-down 
direction of left-going ABS but antiparallel to the spin- 
down direction of right-going ABS. When 9 = 'k j"2., the 
situation is similar. But the phase shift is now the same 
for two ABS because the y-axis makes an angle of 37r/4 
with the spin-down direction of right-going ABS while 
7r/4 with that of left-going ABS. 

When 0 = 7r/4 and 37r/4, the situation is particularly 
interesting. Now the magnetization direction is parallel 
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FIG. 5: (Color online) The CPR of total Josephson current 
for d+p wave with fixed q = 0.56 and various 9: (a) 0, (b) 
7r/4, (c) it/2, (d) 37r/4. In each panel, the results for varying 
values of q from 0 to tt are plotted. Other parameters are the 
same as those in Fig. [S] 

(or antiparallel) to the spin-down direction of one ABS 
but perpendicular to that of the other ABS. So one ABS 
experiences only a phase shift while the other not only 
experiences a phase shift but also shrinks a little. It is 
interesting that the shrinking ABS disappears totally at 
rj = TT, which means that the spin-down particles are 
flipped totally to spin-up and propagate into NCSs as a 
continuous propagating state. For example, when 6 = 
7r/4 and q = tt, the supercurrent is totally carried by 
continuous propagating states instead of discrete ABS for 
the phase difference range tt < ip < 2tt. It is shown that 
the ABS of 0 = 7r/4 and that of 6* = 37r/4 are symmetric 
to each other with respect to the E = 0 axis for the 
reason of symmetry. Similarly, the ABS for (j) = —tt/A is 
symmetric to that for (p = tt/ A with respect to the E = 0 
axis. 

The corresponding total Josephson currents are shown 
in Fig. [5] When 6 = 0, the CPR exhibits a 0-tt tran¬ 
sition with increasing q. The anomalous Josephson ef¬ 
fect does not occur because the two ABS have opposite 
phase shifts as discussed above. When 6 = 7r/4, 9 = tt/2, 
and 6 = 37r/4, the anomalous Josephson current appears. 
The CPR is the same for 9 = tt/A and 9 = "Att/A as the 
symmetry between their ABS shows. And the CPR for 
0 = 7r/2 is only a bit different from that for 9 = tt/A. 
As a whole, the ground-state phase difference is continu¬ 
ously tunable by adjusting q. When 9 = tt/A and q = tt, 
the contribution from cp = tt/A to the supercurrent (not 
shown here) show that the propagating state carried su¬ 
percurrent (0 < V? < TT, see Fig. |4]) equals to that carried 
by the ABS {tt < (p < 2tt). 

Now we discuss the general case of arbitrary q where 
the two gaps are both open. We choose 9 = 7r/4 to let the 
exchange field align with the spin polarization axis of left¬ 
going (right-going) ABS oi <p = tt/A {—tt/A) component. 


FIG. 6: (Color online) The GPR for d-l-p wave with fixed 
6 = tt/A and various q: (a) 7r/4, (b) tt/2, (c) 3tt/A, (d) tt. 
In each panel, q varies from 0 to 1 with a step of 0.2. Other 
parameters are the same as those in Fig. [3] 


The CPRs for various q and q are shown in Fig. [51 The 
case of 77 = 0 is similar to Fig. [3] (a) and thus not shown 
here. For q = 7 r /4 and q = tt/2, the situation is similar 
to that in the case of s-l-p wave. Especially, the case of 
q = tt/ 2 is still important to determine the triplet-singlet 
ratio of NCS because the ground-state phase difference 
ipo varies monotonically with increasing q while the criti¬ 
cal current changes little. What is different from the case 
of s-l-p wave is that ipo is not zero even at g = 0 or 1. That 
is due to the spin-splitting of the momentum factor of gap 
function /(k). For q = Ott/ A and q = tt, the evolution of 
critical current is similar to that in the s-l-p wave case. 
But the evolution of po is very different from that for s-l-p 
wave. Here po changes monotonically with increasing 
q. As discussed previously, that is because the exchange 
field is perpendicular to the spin direction of right-going 
(left-going) ABS oi cp = tt/A {—tt/A) component. These 
ABS with perpendicular spin direction makes difference 
between singlet and triplet pairing. Thus po as well as 
the critical current are different for singlet-dominant and 
triplet-dominant pairing. 


C. d-\-f wave 

For the d+f wave case, the momentum dependence of 
the gap function /(k) = {k^ — ky)/kp. The gap is maxi¬ 
mum at 0 = 0 or tLtt/2. However, the contributions from 
(p = ±7r/2 to the supercurrent are small in comparison 
with that from (p = 0. The contribution from normal 
incidence ky = 0 is still dominant as in the case of s-l-p 
wave. Therefore the situation for d-l-f wave is similar to 
that for s-l-p wave and then the numerical results for d-l-f 
wave are not presented here. 
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V. CONCLUSION 

In summary, we predicted the appearance of anoma¬ 
lous Josephson effect with nonzero ground-state phase 
difference in a NCS/F/NCS junction. The ground-state 
phase difference tpo is proposed to serve as a tool to deter¬ 
mine the triplet-singlet ratio of NCS. The physics picture 
and analytical results are given on the basis of s+p wave, 
while the numerical results and discussion are given on 
both s+p and d+p waves. For s+p wave, ipo reaches 
the extremum when the singlet and triplet components 
have equal magnitude and there is no difference between 
singlet-dominant case and triplet-dominant case gener¬ 
ally. But in the special case of 77 = 7r/2, ipo changes 
monotonically with increasing singlet percentage q. For 
d-fp wave, the monotonic change of ipo with increasing q 


is much more general if only ry is not too small. Interest¬ 
ingly, in the case of d-fp wave, we also find novel states 
in which the supercurrents are totally carried by continu¬ 
ous propagating states instead of discrete ABS. Instead of 
carrying supercurrent, the ABS which here only appear 
above the Fermi energy block the supercurrent flowing 
along the opposite direction. These novel states advance 
the understanding of the relation between ABS and the 
Josephson current. 
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